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摘 要 : 以 2000 年 .2005 年 .2010 年 和 2015 年 ETM 遥感 影像 为 数据 源 ,基于 景观 格局 .生态 敏感 性 和 植被 生产 力 


构建 景观 生态 风险 综合 指数 ,以 ArcGIS 10.2 和 GS +9. 


明 :榆林 市 景观 生态 风险 程度 总 体 较 高 ;2000 年 .2005 年 和 2010 年 以 中 度 风 险 为 主 ,2015 年 以 高 度 风险 为 主 ; 


0 软件 为 平台 ,研究 榆林 市 景观 生态 风险 时 空 分 异 。 结 果 表 


2000—2010 年 ,中 度 风 险 区 范围 在 扩大 ,高 度 风 险 区 范 


甫 在 缩小 ;2010 一 2015 年 ,中 度 风 险 区 范围 迅速 缩小 ,而 高 度 


较 大 。 


风险 区 范围 迅速 扩大 。 榆 林 市 景观 生态 风险 空间 变异 较 小 ,景观 生态 风险 空间 相关 性 中 等 。2000 年 自然 因素 对 景 
观 生 态 风 险 空 间 变异 的 影响 略 大 于 人 为 因素 的 影响 ,2005 年 和 2010 年 的 影响 基本 相当 ,2015 年 则 受 自然 因素 影响 
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生态 风险 评价 是 资源 环境 管理 以 及 生态 建设 的 
重要 依据 ,是 当前 宏观 生态 学 研究 的 热点 之 一 。 它 
起 源 于 20 世纪 80 年 代 , 随 着 生态 风险 评价 研究 的 
快速 发 展 ,评价 尺度 从 生态 系统 向 景观 .区 域 不 断 扩 
大 ,评价 对 象 从 单一 类 型 生态 系统 到 多 种 生态 系统 
的 空间 镶嵌 。 景 观 生态 风 险 指 自然 或 人 为 因素 影响 
下 景观 格局 与 生态 过 程 相互 作用 可 能 产生 的 不 利 后 
FRO 。 其 风险 源 主要 在 于 人 为 景观 的 增加 ,自然 景 
观 的 减少 ,人 类 活动 改变 地 表 生 态 系统 的 格局 和 组 
成 结构 ,从 而 导致 生态 系统 失衡 ,生态 风险 增加 中。 
我 国学 者 对 景观 生态 风险 评价 进行 了 卓有成效 的 研 
究 , 目 前 已 初步 形成 了 具有 一 定 国际 引领 意义 的 景 
观 生态 风险 评价 研究 框架 "。 研 究 案 例 也 比较 多 ， 
研究 对 象 主要 以 流域 “行政 区 “~ 和 城市 地 
域 " ”为 主 , 男 有 海岸 带 及 沿海 地 区 、 海 
FOS Ly Say BE BROT?) ah OP? Ak A) AK AL DK 
WEO 的 景观 生态 风险 评价 研究 ,也 有 针对 具体 建设 
项 目 如 水 电站 建设 的 景观 生态 风险 研究 50 。 

景观 生态 风险 评价 多 以 风险 源 汇 模 式 和 景观 生 
DARIE? I ,尤其 是 景观 生态 模式 应 用 较 广 。 
PUREE A A RWC? 。 景 观 生 态 模式 多 以 景观 指 
数 法 来 刻画 景观 生态 风险 ,所 用 景观 指数 主要 是 反 
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景观 结构 的 指数 ,如 破碎 化 、 分 离 度 指数 .优势 度 
指数 等 。 除 这 两 种 方法 外 ,还 有 用 状态 一 压力 一 
响应 模型 ,结合 投影 寻 踪 法 "" .基于 “潜力 一 连通 
度 一 恢复 力 ” 三 维 准则 建立 评价 指标 体系 、 基 于 
生态 系统 服务 等 评价 景观 生态 风险 。 

目前 景观 生态 模式 着 重 于 从 景观 格局 方面 建立 
评价 指数 ,没有 考虑 景观 生态 系统 的 植被 生产 力 , 植 
被 是 景观 生态 系统 活力 及 其 功能 的 体现 ,植被 生产 
力 高 ,其 景观 生态 系统 活力 强 ,受到 外 界 干 扰 后 , 恢 
复 能 力 强 ,生态 风险 相对 较 小 。 本 文 试图 把 景观 格 
局 和 植被 生产 力 结合 起 来 ,构建 景观 生态 风险 指数 。 

榆林 作为 能 源 化 工 基 地 ,同时 是 陕西 退耕 还 林 
还 草 工 程 的 重点 地 区 之 一 ,人 类 活动 影响 强烈 。 本 
XFJH 2000 年 .2005 年 .2010 年 和 2015 年 4 期 
ETM 遥感 影像 作为 数据 源 ,对 榆林 市 景观 生态 风险 
进行 分 析 评 价 , 以 揭示 研究 区 景观 生态 风险 时 空 变 
化 特征 ,为 该 区 域 生态 环境 管理 和 建设 提供 依据 。 


1 研究 地 区 与 研究 方法 


1.1 研究 区 概况 
榆林 市 (36°57’ ~39°35'N ,107°15’ ~111°15'E) 
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图 1 榆林 两 大 地 貌 单元 


Fig.1 Two main geomorphic units in Yulin 


位 于 陕西 省 最 北部 ,总 面积 43 578 km , 辖 1 区 11 
县 。 地 貌 上 可 分 为 两 大 单元 (图 1) :风沙 滩地 区 和 
黄土 夺 陵 沟 密 区 ,二 者 大 体 以 长 城 为 界 ,毛乌素 沙漠 
南 缘 的 风沙 滩地 区 属 鄂 尔 多 斯 高 原 的 一 部 分 ,黄土 
拷 陵 沟 密 区 属 黄 土 高 原 的 一 部 分 。 气 候 属 温 带 半 干 


107°E 108°E 109°E 110°E 111°E 
le T T T T 


39°N 


38°N 


37°N 


107°E 108°E 109°E 110°E 111°E 


39°N 


ony) Mil 中 覆盖 度 章 地 
38 N | mmm 高 覆盖 度 草地 


37°NE 是 eee 
0 30 60 90 120km (c) 2010 年 


旱 大 陆 性 季风 气候 ,四 季 分 明 ,年 平均 气温 10 C ,年 
平均 降水 量 400 mm 左右 ,年 内 降水 量 多 集中 于 6 一 
9 月 。 土壤 主 要 有 风沙 土 、 黄 绵 土 .新 积 土 . 灰 钙 土 、 
REE RIE . 盐 圭 等。 地带 型 植被 为 干草 原 和 和 森 
林 草 原 , 主要 植被 有 和 森林、 干草 原 EAA, BE fi] VD AE 
植被 . 盐 生 植被 等 。 榆 林 矿 产 资源 极其 丰富 ,有 煤 、 
天 然 气 石油 . 盐 等 ,为 陕 北 能 源 化 工 基地 的 核心 区 
域 。 

1.2 资料 与 方法 

1.2.1 资料 与 数据 处 理 ”本文 所 用 遥感 影像 资料 
SURF FREE SHE z (http://www. csdb. en), 在 
2000 年 .2005 年 .2010 年 和 2015 年 ETM 遥感 影像 
解 译 得 到 的 4 期 景观 类 型 (景观 类 型 是 参考 中 国 科 
学 院 土地 覆盖 分 类 体系 ,结合 研究 需要 划分 ) 图 ( 解 
译 精度 均 在 87.5% 以 上 ,图 2) 的 基础 上 ,从 两 个 方 
面 计算 景观 生态 风险 指数 ,一 是 以 榆林 市 为 整体 景 
观 ,计算 4 个 时 相 的 景观 生态 风险 指数 ,以 确定 整个 
榆林 景观 生态 风险 大 小 和 动态 变化 ;二 是 将 研究 区 
域 用 10 km x 10 km 的 网 格 进行 划分 , 共 分 为 527 块 
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2 榆林 市 各 年 份 景观 生态 


Fig.2 Landscape ecology in Yulin in different years 


THR 


景观 生态 小 区 ,计算 4 期 每 块 生态 小 区 的 景观 生态 
风险 指数 。 

地 统计 学 方法 是 分 析 区 域 化 变量 空间 异 质 性 的 
有 效 工具 , 它 是 通过 变异 隐 数 确定 变量 因子 的 空间 


研 R 


算 之 后 ,再 将 其 用 公式 (3 ) 标 准 化 ,以 消除 其 量 纲 ， 
之 后 用 公式 (1) 来 计算 景观 生态 风险 指数 。 
X* = (X -X nin) / (Xmax — Xin) (3) 


7S Sp AR BE DY 25 BRR BE OY ,变异 函数 是 地 统计 学 
特有 的 基本 工具 。 本 文 利 用 SPSS 20.0 对 各 年 景观 
生态 风险 小 区 的 风险 指数 进行 描述 性 统计 和 正 态 分 
布 检验 ,采用 GS + 9.0 软件 进行 地 统计 分 析 , 建 立 
变异 函数 的 半 方 差 理 论 模型 ,研究 其 空间 分 异 。 
1.2.2 景观 生态 风险 评价 ”景观 生态 系统 结构 越 
复杂 ,多样 性 越 高 , 景观 类 型 越 丰富 ,抵御 风险 的 能 
力 越 强 ”2 。 景 观 破碎 化 是 景观 生态 系统 生物 多 样 
性 丧失 的 一 个 主要 原因 ,因此 ,景观 破碎 化 也 会 加 大 
景观 生态 风险 。 景 观 连 接 度 高 ,景观 要 素 间 的 相互 
联系 多 ,其 风险 小 。 敏 感性 高 的 景观 ,其 生态 系统 容 
易 受 损 ,其 风险 也 高 。 植 被 是 陆地 景观 生态 系统 活 
力 的 象征 ,景观 生态 系统 植被 生物 量 大 、 履 盖 度 高 ， 
恢复 能 力 强 ,生态 风险 小 。 因 此 ,本 文 从 景观 破碎 
化 .景观 连接 度 .景观 生态 敏感 性 .景观 多 样 性 以 及 
植被 生产 力 角 度 来 刻画 景观 生态 风险 。 本 文采 用 景 
观 破碎 化 指数 (FRAG) iE Be BE FS BX (CONNECT) 、 
生态 敏感 性 指数 (SENS1) OWL FE VETS BC ( SHDI) 
和 NDVI 来 定量 表达 景观 生态 风险 ,其 公式 为 : 

ILER =1 - aCONNECT +BFRAG - pSHDI - 
SNDVI + eSENSI, (O<ILER<1) (1) 


AP LER 为 景观 生态 风险 指数 ;a、B6、p.6 和 分 
别 为 相应 项 的 权重 ,用 IDRISI 软件 中 的 权重 计算 模 
块 通过 两 两 比较 ,构建 判断 矩阵 来 计算 ,其 值 分 别 为 
0.221 6 .0.221 6.0.2216.0.2716 和 0.063 6。 

CONNECT 和 SHDI 为 景观 格局 指数 软件 
FRAGSTATS 中 景观 水 平 上 的 两 个 指数 。 

NDVI 利用 250 m 2000 年 .2005 年 .2010 年 和 
2015 年 年 平均 MODIS NDVI 值 计 算 ; 

笔者 提出 一 个 新 的 景观 破碎 化 指数 FRAG, H 
计算 公式 如 下 : 

FRAG= (NPxTE)/(TAx10000) (2) 


式 中 ;NP 是 景观 内 斑 块 总 数 ; TE 为 景观 内 边界 总 
长 ;74 为 景观 总 面积 。 

景观 生态 敏感 性 指数 SENSI 从 3 个 方面 来 计 
S 

人 性。 具体 计算 方法 见 文献 L35] 。 

CONNET PR4C、SHDTNDJ 和 SENSI 参数 计 


SF 


式 中 : 为 参数 标准 化 后 的 值 ,和 是 每 个 样本 (每 个 
生态 小 区 ) 的 相应 参数 值 ,X,,. 和 XX,;, 分 别 是 所 有 样 
本 中 相应 参数 的 最 大 值 和 最 小 值 。 

景观 生态 风险 指数 计算 完成 并 对 其 进行 地 统计 
检验 后 ,采用 ArcGIS 9. 3 地 统计 模块 中 普通 Krige 
插值 法 进行 空间 插值 ,按照 景观 生态 风险 指数 的 大 
小 及 分 布 情况 ,把 景观 生态 风险 分 为 4 级 : 低 度 风险 
(ILER <0.4) PERK (0. 4<ILER <0. 6) ,高度 风 
K (0. 6 <ILER <0.8) 和 极 高 度 风险 (ILER 二 0. 8) , 
以 说 明 景 观 生态 风险 的 相对 大 小 ,得 到 榆林 市 景观 
生态 风险 分 级 图 ,并 研究 景观 生态 风险 级 别 的 时 空 
变化 。 

景观 生态 风险 各 等 级 之 间 的 转化 用 IDIRSI 软 
件 中 的 CROSSTAB 模块 进行 运算 ,生成 转化 图 , 统 
计 各 等 级 之 间 的 转化 面积 ,人 然后 计算 转 出 率 和 转 入 
率 。i 级 景观 生态 风险 区 转化 为 j 级 景观 生态 风险 
区 的 转 出 率 为 : 


ROUT ea 100% , (iF; 4 
taj AREA, x o, (iAj) (4) 
级 景观 生态 风险 的 转 出 率 为 : 

ROUT, = Ý, ROUT,,, (ij) (5) 


由 i 级 景观 生态 风险 区 转 入 到 j 级 景观 生态 风 


险 区 的 转 和 人 率 为 : 
AREA,,» wae 
RIN = AREAS x 100% , (i#j) (6) 
级 景观 生态 风险 区 的 转 入 率 为 : 
RIN, = È, RING, (i #3) (7) 


(4) ~ (7) sth, ROUT, ,为 某 一 时 期 i 级 景观 生态 
风险 区 转化 为 7 级 景观 生态 风险 区 的 转 出 率 ;ARE- 


hj 为 某 一 时 期 i 级 景观 生态 风险 区 转化 为 7 级 景 


观 生态 风险 区 的 面积 (km ) ;4RE4, 为 某 一 时 期 期 
初 i 级 景观 生态 风险 区 的 面积 (km? ) ; ROUT, 为 某 
一 时 期 i 级 景观 生态 风险 区 的 转 出 率 ; RING ) 为 某 
一 时 期 i 级 景观 生态 风险 转 为 j 级 景观 生态 风险 区 
的 转 信 率 ;ARE4, 为 某 一 时 期 期 末 j 级 景观 生态 风 
险 区 的 面积 (km ) ; RIN, 为 某 一 时 期 j 级 景观 生态 
风险 区 的 转 入 率 。 
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2 结果 与 讨论 


2.1 地 统计 分 析 

2.1.1 景观 生态 风险 指数 统计 特征 表 1 是 4 期 
527 块 景观 生态 小 区 景观 风险 指数 的 统计 特征 值 。 
可 以 看 出 ,2000 年 2005 年 和 2010 年 的 均值 都 接近 
0.5 ,2015 年 均值 稍 高 ,为 0.685。 从 表 1 中 的 偏 度 
和 峰 度 系数 可 见 研 究 区 景观 生态 风险 指数 呈正 态 
分 布 。 


表 1 研究 区 景观 风险 指数 统计 特征 值 
Tab.1 Statistical characteristic values of landscape 


risk index in the study area 


年 份 ” 中 位 数 ”均值 标准 差 偏 度 GE ”方差 


2000 0.500 0.501 0.118 一 0. 364 1.607 0.014 
2005 0. 500 0.501 0.118 一 0. 364 1.607 0.014 
2010 0. 464 0.476 0. 103 0. 884 5.376 0.011 
2015 0. 691 0.685 0. 142 —0. 682 1.342 0.020 
2.1.2 景观 生态 风险 空间 变异 结构 特征 ”本 研究 


的 一 半 , 且 各 年 景观 生态 风险 样本 数据 (527 个 生态 
风险 小 区 的 景观 生态 风险 指数 ) 呈正 态 分 布 ,满足 
地 统计 分 析 要 求 。 

从 2000—2015 年 ,景观 生态 风险 的 块 金 值 都 比 
较 小 ( 表 2 ,图 3) ,说 明 随机 部 分 的 景观 生态 风险 空 
间 蜡 质 性 比较 小 ,而 自 相关 变异 较 大 。 基 台 值 也 都 
比较 小 ,表明 榆林 市 景观 生态 风险 的 空间 异 质 性 总 
体 来 讲 比较 低 。 相 对 来 说 ,2015 年 景观 生态 风险 总 
的 空间 异 质 性 程度 最 高 ,而 2010 年 总 的 空间 异 质 性 
最 低 。2000 年 .2005 年 .2010 年 和 2015 年 景观 生 
oil cn clone ea ad 

分 别 占 总 空间 异 质 性 的 41. 240% 、49. 752% 、 
49.574% 和 34.526% , 且 主 要 体现 在 10 km UFR 
度 上 ;而 由 空间 自 相 关 引 起 的 空间 蜡 质 性 分 别 占 总 

空间 异 质 性 的 58. 760% 、50. 248% 、50. 426% 和 
65.474% ,分 别 体现 在 10.0 ~57.4 km 10.0 ~ 835. 2 
km 10. 0 ~ 294.3 km 和 10.0 ~ 122.4 km 范围 内 。 
2005 年 和 2010 年 由 随机 部 分 引起 的 空间 异 质 性 程 


中 步 长 设 为 10 km, $k 22 组 ,小 于 样 点 间 最 大 距离 


表 2 榆林 市 景观 生态 风险 半 变 异 函 数 


度 比 2000 年 高 。 可 见 ,2000 年 资源 开发 .退耕 还 林 


理论 模型 参数 


Tab.2 Theoretical model of semivariograms and parateters of landscape ecological risk in Yulin 


年 份 模型 块 金 值 基 台 值 块 金 值 / 基 台 值 /% 变 程 /km r 残 差 ( x106) 
2000 球状 模型 0.005 36 0.013 32 41.240 57.4 0.874 9.211 
2005 指数 模型 0.010 03 0. 020 16 49.752 835.2 0.764 16.630 
2010 站 数 模型 0.005 82 0.011 74 49.574 294.3 0. 890 5.391 
2015 球状 模型 0.007 14 0. 020 68 34. 526 122.4 0.979 7.570 
0.0147, (a) 2000 年 eee 0.0156, (Œ) 2005 年 ee 
0.011 0 0.0117 ae 
WW rm ees 
$R 0.0073 IR 0.0078} o 
at at i 
0.003 7 | 0.003 9 + 
0 112 000 224 000 0 112 000 224 000 
0.0112 (c) 2010 年 0.0215 (d) 2015 年 
0.008 4 0.016 1 
x ad 
R 0.0056 F, 要 0.0107 
+ [ m 
0.002 8 上 0.005 4 + 
0 112 000 224 000 0 112 000 224 000 
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图 3 榆林 市 不 同年 份 景 观 生态 
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Fig.3 Semivariograms for landscape ecological risk in Yulin of different years 
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等 人 为 因素 对 景观 生态 风险 的 空间 变异 影响 略 小 于 
由 地 形 .气候 等 自然 因素 的 影响 ,而 2005 年 和 2010 
年 ,人 为 因素 和 自然 因素 的 影响 基本 相当 ,2015 年 
空间 异 质 性 主要 由 自然 因素 引起 。 

See A (Co) 与 (Co +C) 基 台 值 之 比 可 以 表明 区 
域 化 变量 空间 相关 性 程度 , 当 Cy (Cy +C) < 25% 
时 ,变量 具有 强烈 的 空间 相关 性 ;在 25% ~75% 时 ， 
区 域 化 变量 具有 中 等 的 空间 相关 性 ; 当 C6/(Co+ 
C) >75% 时 ,变量 空间 相关 性 很 弱 ”” 。 从 块 金 值 
与 基 台 值 之 比 来 看 ,这 4 a 榆林 市 景观 生态 风险 的 
块 金 值 和 基 台 值 之 比 介 于 25% ~75% 之 间 , 因 此 ， 
景观 生态 风险 空间 相关 性 中 等 。 
2.2 景观 生态 风险 时 间 变 化 

就 整个 榆林 市 而 言 ,2000 年 .2005 年 .2010 年 
和 2015 年 的 景观 生态 风险 指数 分 别 为 0. 752、 
0.713 .0.685 和 0. 792, K 2010 年 外 ,景观 生态 风 
念 指数 均 大 于 0.7 ,说 明 榆 林 市 景观 生态 风险 比较 
大 ;从 时 间 变 化 上 来 看 ,2000 一 2010 年 景观 生态 风 
险 指 数 在 下 降 ,2010 一 2015 年 景观 生态 风险 指数 在 
增 大 。 

从 风险 等 级 的 时 间 动 态 变化 来 看 ( 表 3) ,2000 
年 .2005 年 和 2010 年 ,榆林 以 中 度 风险 为 主 ,其 面 
积 均 占 到 整个 榆林 市 的 70% 以 上 ;2015 年 以 高 度 风 
险 为 主 ; 2005 年 和 2010 年 没有 极 高 度 风 险 区 。 
2010 一 2015 年 ,中 度 风 险 区 面积 迅速 减 小 , 而 高 度 
风险 区 面积 迅速 增加 ,同时 极 高 度 风 险 区 面积 也 在 


BS 


究 


研 


RI 榆林 各 级 景观 生态 风险 占 全 区 面积 比 


Tab.3 The area proportions of landscape ecological risk 


at all levels in the whole area in Yulin 1% 
风险 级 别 ” 等 级 ”2000 年 ”2005 年 ”2010 年 ”2015 年 
低 度 风险 1 0 1.438 1.049 0 
中 度 风险 2 79.066 92.614 97.951 9.536 
高 度 风 险 3 20. 882 5.948 1.000 80. 046 
极 高 度 风 险 4 0.052 0 0 10. 418 
合计 100 100 100 100 


积 有 所 增加 ,高 度 风 险 区 面积 有 所 减 小 。2000 一 
2015 年 极 高 度 风险 区 面积 扩大 近 10.4% ,高 度 风险 
区 面积 扩大 近 60% ,中 度 风险 区 面积 缩小 较 大 。 

从 各 级 别 的 风险 转化 情况 来 看 ( 表 4 ) ,2000 一 
2005 年 , 极 高 度 风险 区 全 部 转 为 低 度 风 险 区 ,其 次 
是 高 度 风险 区 , 转 出 率 为 80. 414% ,主要 转向 了 中 
度 风险 区 ; 低 度 风险 区 转 入 率 最 大 ,为 100% ,主要 
由 中 度 和 高 度 风 险 区 转化 而 来 。2005 年 与 2000 年 
相 比 , 极 高 度 风险 区 完全 消失 (图 4)。 从 2005 一 
2010 年 , 低 度 风险 区 转 出 率 高 ,为 100% ,转向 了 中 
度 风 险 区 ,其 空间 位 置 完全 发 生 了 变化 。 其 次 为 高 
度 风险 区 ,主要 转向 了 中 度 风险 区 ; 低 度 风险 区 的 转 
入 率 最 高 ,接近 100% ,主要 由 中 度 风险 区 转化 而 
来 ,其 次 为 高 度 风险 区 ( 转 入 率 为 38.078%)。 从 
2010—2015 年 , 低 度 风险 区 全 部 转 出 ,主要 转 为 中 
度 风险 区 ,其 次 是 中 度 风 险 区 , 转 出 率 为 91.512% ， 
主要 转向 了 高 度 风 险 区 ; 转 和 人 率 以 极 高 度 风 险 区 和 


增加 。2010 年 比 2000 年 和 2005 年 中 度 风 险 区 面 


> 


高 度 风险 区 最 高 ,分 别 达到 100% 和 98. 966% ,主要 


表 4 2000 一 2010 年 景观 生态 风险 各 等 级 之 间 的 转化 


Tab.4 Transition ratios of different landscape ecological risk levels in Yulin from 2000 to 2010 /% 
转换 类 型 2000 一 2005 年 2005 一 2010 年 2010 一 2015 年 2000 一 2015 年 
转 出 率 RAR 转 出 率 转 入 率 转 出 率 转 人 率 转 出 率 RAR 

2 一 ] 1.406 77. 324 1.131 99.855 0 0 0 0 

31 1.561 22.676 0 0 0 0 0 0 

1—2 = = 100. 000 1.466 84.115 9.256 = = 

3 一 2 78.853 17.779 89. 588 5.440 0. 346 3.553 1.775 3.887 

4—2 100.000 0.057 = = = = 0 0 

43 0 0 = z z T 100. 000 0. 065 

2 一 3 2.350 31.237 0.411 38.078 80.876 98.966 81.637 80. 637 

3—4 0 0 0 0 0 0 24. 252 48. 609 

2—4 0 0 0 0 10.636 100.000 6.772 51.391 

1 一 3 = = 0 0 15. 885 0.208 - - 
注 :2 一 1 为 中 度 风险 区 转化 为 低 度 风险 区 ;3 一 1 为 高 度 风险 区 转化 为 低 度 风险 区 ;1 一 2 为 低 度 风险 区 转化 为 中 度 风险 区 ;3 一 2 为 高 度 风险 区 
转化 为 中 度 风险 区 ;4 一 2 为 极 高 度 风险 区 转化 为 中 度 风险 区 ;4 一 3 为 极 高 度 风险 区 转化 为 高 度 风险 区 ;2 一 3 为 中 度 风险 区 转化 为 高 度 风险 
区 ;3 一 4 为 高 度 风险 区 转化 为 极 高 度 风 险 区 ;2 一 4 为 中 度 风险 区 转化 为 极 高 度 风 险 区 ;1 一 3 为 低 度 风险 区 转化 为 高 度 风险 区 ，。 
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4 榆林 市 不 同年 份 景观 生态 风险 分 布 


Fig.4 Distribution of landscape ecological risk in Yulin in different years 


由 中 度 风 险 区 转化 而 来 。2000 一 2015 年 , 极 高 度 风 
险 区 全 部 转 为 高 度 风 险 区 ,81. 637% 的 中 度 风 险 区 
转 为 高 度 风险 区 ; 极 高 度 风 险 区 的 转 人 率 也 达到 
100% ,其 中 48. 609% 由 高 度 风 险 区 转化 而 来 ， 
51. 391% 由 中 度 风 险 区 转化 而 来 。 除 此 而 外 ,高 度 
风险 区 的 转 入 也 很 高 ,达到 80.702% 。 

榆林 市 景观 生态 风险 及 其 变化 受 自然 条 件 及 人 
类 活动 的 共同 影响 。 自 然 方面 主要 受气 候 、 地 形 A 
被 及 地 表 组 成 物质 等 的 影响 。 榆 林 市 地 处 半 干 旱 的 
温带 ,降水 稀少 ,水 资源 短缺 ,西北 部 为 风沙 草滩 区 ， 
东南 部 为 黄土 丘陵 沟 密 区 ,生态 脆弱 ,加 之 区 内 矿产 
资源 丰富 ,资源 开发 等 人 为 活动 强烈 ,因而 ,其 整体 
景观 生态 风险 较 高 。 榆 林 气温 呈现 增加 趋势 ”3， 
而 降水 呈现 不 稳定 性 ,这 会 导致 景观 生态 风险 的 
增强 。 


人 为 活动 方面 , 自 2003 年 原 国家 计 委 批准 陕 北 
能 源 化 工 基 地 在 榆林 启动 建设 以 来 ,榆林 市 发 展 加 
速 , 已 经 形成 榆林 和 神 府 2 个 经 济 开发 区 , 榆 神 煤 化 
工 、 榆 横 煤 化 工 等 6 个 工业 园区 。2000 一 2010 年 ， 
累计 建成 各 类 重点 项 目 70 个 ,完成 投资 1.26 x 10" 
元 。2010 年 以 后 能 源 化 工 基地 继续 建设 ,2010 一 
2015 年 ,建成 转化 产能 煤 基 甲醇 2.20 x 10° t, HER 
甲醇 醋酸 2. 00 x 10° t、 煤 制 油 品 1.60 x 10° t、 聚 氧 
乙烯 1.35 x 10° t R ZG. 20 x105t 聚 丙烯 1.20 x 
10° t 兰 炭 5.00 x10”t。 这 些 工业 项 目 大 都 以 煤化 
工 \ 盐 化 工 、 煤 电 等 高 耗 水 行业 为 主 ,而 全 市 水 资源 
总 量 约 32 x10 m ,人 均 水 资源 865 mi ,和 远 低 于 陕西 
省 和 全 国平 均 水 平 , 属 重度 缺 水 地 区 ,加 之 原煤 、 原 
油 等 工业 品 产量 逐年 增长 ( 表 5) ,因而 更 加 剧 了 水 
资源 的 消耗 ,使 地 表 水 减少 ,地 下 水 系 破 坏 , 地 下 水 


TF 


位 下 降 ,湿地 面积 萎缩 ,从 而 增 大 景观 生态 风险 。 如 
仅 神 木 县 境内 就 有 10 条 河 断 流 .20 多 眼 泉 井 干 枯 ， 
乌 兰 木 伦 河 的 河床 被 占 去 90% 。 全 榆林 市 湖泊 由 
开发 前 的 869 个 减少 为 79 个 。2000 年 能 源 化 工 基 
地 建设 尚未 启动 ,因此 ,景观 生态 风险 的 空间 异 质 性 
主要 由 自然 因素 引起 。 而 到 2005 年 和 2010 年 , 随 
着 能 源 化 工 基地 建设 启动 和 实施 ,许多 工矿 项 目 入 
驻 , 人 为 影响 逐渐 增 大 ,随机 因素 引起 的 空间 变异 接 
近 50% ,这 时 人 为 因素 和 自然 因素 对 空间 异 质 性 的 
影响 相当 。 到 了 2015 年 ,由 随机 因素 引起 的 空间 异 
质 性 下 降 到 34.526% ,可 能 是 工农 业 生 产 活动 的 增 
强 导 致 了 某 些 自然 因素 ,如 水 ,成 为 影响 景观 生态 风 
险 空 间 变 异 的 主要 因素 ,关于 这 一 点 有 待 进一步 深 
入 研究 。 

退耕 还 林 还 草 也 是 人 为 活动 的 一 个 主要 方面 。 
1999—2006 年 是 榆林 退耕 还 林 还 草 工程 实施 的 主 
要 时 期 ,2006 年 之 后 ,退耕 还 林 还 草 规 模 逐 渐 减 小 ， 
到 后 期 耕地 数量 不 降 反 增 ( 表 6) ,林地 在 2005 年 后 
逐渐 减少 ,牧草 地 在 2000 一 2013 年 期 间 一 直 在 减 
少 , 而 建设 用 地 和 工矿 用 地 一 直 在 增加 ,在 短 短 的 
2010 一 2013 年 的 3 a 中 ,二 者 的 增加 量 比 2000 一 
2005 年 间 的 增加 量 还 要 大 ,后 期 耕地 、 建 设 用 地 及 
工矿 用 地 的 增加 无 疑 会 增加 景观 破碎 化 ,导致 景观 
生态 敏感 性 增强 。 初 期 ,退耕 还 林 还 草 工程 实施 ， 
工业 园区 建设 等 活动 处 于 建设 期 ,引进 的 项 目 不 多 ， 
建设 活动 相对 来 说 比较 少 , 加 之 好 多 项 目 还 没 投产 ， 
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风险 逐年 减 小 ;但 后 期 随 着 退耕 还 林 还 草 工 程 实施 
减缓 ,或 任务 完成 ,入 驻 工 业 园区 的 企业 越 来 越 多 ， 
生产 规模 越 来 越 大 ,产能 越 来 越 大 ,对 水 资源 的 消耗 
也 会 越 来 越 大 ,水 资源 矛盾 越 来 越 突 出 ,如 2004 年 
工业 用 水 量 为 0. 529 x 10° m ,2010 年 为 1.297 x 
10° m° ,到 2013 年 为 1.796 xl0sm ,除了 工业 用 水 
外 ,农业 也 是 用 水 大 户 ,2005 年 后 ,耕地 不 断 增加 ， 
无 疑 增 加 了 水 资源 的 消耗 ,2004 年 榆林 农业 用 水 量 
4.927 x 10ë m° ,2010 年 为 4. 920 x 10° mi ,到 2013 
年 为 5.296 x 10° m ™™ ,生态 环境 的 变化 具有 一 定 
的 滞后 性 ,因此 ,2000 一 2010 年 生态 风险 指数 减 小 ， 
2015 年 生态 风险 指数 又 增 大 。 
2.3 ”景观 生态 风险 空间 分 异 

从 景观 生态 风险 的 空间 差异 来 看 ,2000 年 、 
2005 年 和 2010 年 ,风沙 滩地 区 和 黄土 丘陵 区 都 以 
中 度 风险 为 主 , 风 沙滩 地 区 这 3 a 中 度 风 险 区 占 到 
该 区 面积 的 96% 以 上 ( 表 7) ,2010 年 中 度 风 险 区 面 
积 接近 100% 。 黄 土 丘陵 沟 餐 区 中 度 风 险 面 积 从 
2000 年 的 67. 453% ,到 2005 年 的 88.715% ,再 到 
2010 年 的 96.593% ,逐渐 增 大 ,到 2010 年 已 占 绝对 
优势 。2015 年 ,空间 分 布 格局 出 现 了 明显 变化 , 黄 
土 丘 陵 区 和 风沙 滩地 区 均 以 高 度 风 险 为 主 ,不 仅 如 
此 ,2 个 区 域 均 出 现 了 10% 以 上 的 极 高 度 风险 区 。 

风沙 滩地 区 气候 为 温带 半 干 旱 气 候 , 年 降水 量 
多 在 400 mm 以 下 , 蒸发 量 在 1 000 mm 以 上 , 地 貌 


表 5 榆林 部 分 工业 品 产 量 
Tab.5 Outputs of some industrial products inYulin 
原煤 原油 原盐 洗煤 精 甲 醇 RAC Mi “RR 天 然 气 发 电量 
年 份 8 m? 
/10 t /108 mi /108 kW h 
2000 3 125.0 89.00 7.50 一 = 一 = 34. 00 = 
2005 6 560.3 476.40 26.30 = 46. 88 = = 60. 50 68.45 
2010 25 732.0 983.00 40.98 = 119. 86 16.41 959.56 109. 86 356. 96 
2015 36 103.5 1 186.94 134. 80 7 053.05 222.07 131.36 2 499.11 151.13 629. 34 
TE :数据 来 源 于 榆林 市 国民 经 济 和 社会 发 展 统计 公报 。 
6 2000—2013 年 主要 土地 利用 类 型 的 变化 
Tab.6 Land use change in Yulin from 2000 to 2013 /hm? 
时 段 耕地 园地 林地 牧草 地 建设 用 地 [ 矿 用 地 荒草 地 盐碱地 沙 地 
2000—2005 年 -62 093.98 23 804.79 176 115.506 -40976.75 13 687.00 739.37 10 736.94 一 233.40 一 9 092.60 
2005—2010 4Æ 54569.38 93 674.91 -31 493.513 -349 513.86 28723.00 6 098.57 256 490.75 -8 658.47 -33 416.50 
2010—2013 年 5 916. 16 480.35 -8 881.950 -5 623.61 39 686. 13 6611.16 -410 726.81 - 502.54 -2 618.72 


注 : 负 值 表示 减少 ， 


其 余 为 增加 ; 数据 来 源 于 2000 年 .2005 年 .2010 年 和 2013 年 榆林 市 土地 变更 调查 数据 。 
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表 7 黄土 丘陵 沟 直 区 和 风沙 滩地 区 各 级 景观 生态 风险 占 其 面积 比 


Tab.7 The area proportions of landscape ecological risk at all levels in the loess hilly region and 


aeolian sandy region /% 
风险 等 级 2000 年 2005 年 2010 年 2015 年 

工区 I K 工区 [I 区 工区 工区 工区 工区 
1 0 0 2.392 0 1.743 0. 006 0 0 

2 67.455 96.551 88.715 98. 488 96. 593 99.994 15.395 0.712 

3 32.458 3.449 8. 893 1.512 1. 664 0 73.991 89. 165 

4 0.087 0 0 0 0 0 10.614 10. 123 
合计 100 100 100 100 100 100 100 100 


注 ; [区 为 黄土 丘陵 沟 蜜 区 ; I Ky a H 


区 


由 东 向 西 主要 是 沙漠 和 和 干 湿 沙 滩 , 地 带 性 植被 以 温 
带 干 草原 和 荒漠 草原 .沙化 植被 为 主 ,土壤 以 风沙 土 
为 主 , 地 下 水 丰富 ,地 下 水 埋 深 5 ~50 m。 有 丰富 的 
煤炭 天 燃气 和 盐 类 资源 ,近年 来 也 是 耕地 补充 的 主 
要 地 区 之 一 ,区 内 有 靖 边 能 源 化 工 综合 利用 产业 园 
区 \ 榆 横 工 业 区 \ 榆 神 工业 区 等 ,因而 人 为 活动 强烈 ， 
风险 较 高 。 黄 土 丘陵 沟 歼 区 气候 为 暖 温带 半 和 干旱 气 
候 , 降 水 量 400 ~ 500 mm ,降水 主要 集中 在 夏季 , 且 
多 大 十 和 暴雨 ;地 犁 以 操 入 和 沟谷 为 主 ,地 表 破 碎 ， 
起 伏 较 大 ,地 面 坡度 大 ;地 下 水 埋 深 150 ~ 250 m, , 严 
重 缺 水 ,只 有 在 河谷 川 道 , 地 下 水 丰富 , 埋 深 50 ~ 
150 m; 生 态 脆 弱 , 区 内 矿产 主要 有 天 然 气 、 岩 盐 、 煤 
等 ,区 内 有 绥 ( 德 ) 米 ( 脂 ) 佳 (县 ) 盐 化 工业 园区 R 
堡 煤 焦化 工业 园区 等 ,人 为 活动 也 较 强 烈 ,生态 风险 
较 高 。 可 见 ,无 论 是 风沙 滩地 区 和 黄土 丘陵 沟 餐 区 ， 
本 身 自然 条 件 决定 了 其 生态 风险 较 高 ,加 之 ,矿产 资 
源 开发 等 人 为 活动 强烈 ,使 生态 风险 更 高 。2 个 区 
的 生态 风险 时 间 动 态 变化 主要 还 是 人 为 因素 引起 
的 ,根据 2013 年 榆林 气候 公报 数据 (《2013 年 气候 
公报 》, 榆 林 市 气象 局 ,2014 年 1 A) ,从 2000 一 2013 
年 位 于 风沙 滩地 区 的 最 大 胡 泊 一 一 红 碱 淖 水 体面 积 
一 直 旦 减少 趋势 ,即使 在 降水 增多 的 年 份 , 红 碱 淖 水 
体面 积 也 在 减少 ,说 明 红 碱 淖 水 体面 积 减少 不 是 气 
候 因素 ,而 是 人 为 因素 。 
2.4 研究 的 局 限 性 

本 文 只 是 研究 了 4 个 年 份 的 情况 ,如 果 能 逐年 
研究 景观 生态 风险 ,能 更 好 地 反映 其 时 间 动 态 变化 。 
水 是 影响 榆林 市 景观 生态 的 一 个 很 重要 的 因素 ,本 
文 只 是 从 景观 格局 和 NDVI 角度 来 刻画 景观 生态 风 
险 ,虽然 NDVI 是 气候 .水 、 土 等 的 综合 反映 ,但 
如 果 能 够 将 水 文 因素 纳入 到 景观 生态 风险 评价 中 ， 
会 使 研究 更 加 完善 ,这 也 是 今后 努力 的 方向 。 


3 结论 


< 


(1) 从 2000—2010 年 ,景观 生态 风险 指数 在 减 
小 ,而 2010 一 2015 年 景观 生态 风险 指数 在 增 大 。 
2000 年 .2005 年 和 2010 年 ,以 中 度 风险 为 主 ,2015 
年 以 高 度 风 险 为 主 。2000 一 2010 年 中 度 风险 面积 
在 扩大 ,高 度 风 险 区 面积 在 减 小 ,2010 一 2015 年 高 
度 风 险 区 面积 增加 ,中 度 风险 区 面积 在 减 小 。 榆 林 
市 地 处 半 和 干旱 的 温带 ,降水 稀少 ,水 资源 短缺 ,生态 
脆弱 ,人 类 活动 强烈 ,总体 来 说 景观 生态 风险 较 高 。 
退耕 还 林 还 草 工程 的 实施 使 景观 生态 风险 降低 ,能 
源 化 工 基 地 生产 规模 扩大 ,生产 能 力 提高 ,使 景观 生 
态 风险 增 大 。 

(2) 榆林 市 景观 生态 风险 空间 变异 较 小 ,景观 
生态 风险 空间 相关 性 中 等 。 由 随机 因素 引起 的 空间 
分 异 主要 体现 在 10 km 范围 内 ,而 由 结构 性 因素 引 
起 的 分 异 主 要 体现 在 10 km 以 外 范围 。 气 候 、 地 貌 、 
土壤 等 自然 因素 主要 引起 景观 生态 风险 大 范围 的 空 
间 变 异 , 人 为 因素 主要 引起 小 范围 内 的 空间 变异 。 
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Landscape Ecological Risk and Its Spatiotemporal Variation in Yulin 


SHI Yu-qiong'*”, WANG Ning-lian?”, LI Tuan-sheng*, WANG Han’ 
(1. College of Forestry ,Shanxi Agricultural University , Taigu 030800 , Shanxi , China ; 
2. College of Urban and Environmental Sciences ,Xi’an 710124 ,Shaanxi , China; 
3. Shaanxi Province Key Laboratory of Surface System and Environmental Carrying Capacity ,Xi’an 710124 , Shaanxi , China ; 
4. College of Earth and Resources Sciences ,Chang’an University ,Xi’an 710054 , Shaanxi , China) 


Abstract; The ETM remote sensing images in 2000 ,2005 ,2010 and 2015 were taken as the data sources to con- 
duct the comprehensive index of landscape ecological risk in Yulin based on the landscape pattern , ecological sensi- 
tivity and vegetation productivity. The ArcGIS 10.2 and GS +9.0 software were used as the platforms to study the 
spatiotemporal differentiation of landscape ecological risk in the study area. The results showed that the landscape 
ecological risk in Yulin was holistically high. Risk of the main part of the regional landscape ecological ecosystem 
was moderate in 2000 ,2005 and 2010, but high in 2015. From 2000 to 2010, the area with moderate risk was en- 
larged , but that with high risk was reduced. From 2010 to 2015 ,the area with moderate risk was reduced rapidly , 
but that with high risk was enlarged rapidly. The spatial variation of landscape ecological risk was low in Yulin, and 
the spatial correlation of landscape ecological risk was moderate. The effect of natural factors (such as landform and 
climate ) on the spatial variation of landscape ecological risk was slightly higher than that of human factors in 2000 , 
and it was almost equal to each other in 2005 and 2010. In 2015 , however, the effect of natural factors on the spatial 
variation of landscape ecological risk was high. 

Key words; landscape ecological; risk index; spatiotemporal variation; structural feature; Yulin; Shaanxi Prov- 


ince 


